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Abstract 

Background: In Zambia, the presence of bovine tuberculosis in both wild and domestic animals has long been 
acknowledged and mutual transmission between them has been predicted without any direct evidence. 
Elucidation of the circulating Mycobacterium bovis strains at wild and domestic animals interphase area in Zambia, 
where bovine tuberculosis was diagnosed in wildlife seemed to be important. 

Results: A PCR identified 15 and 37 M. bovis isolates from lechwe and cattle, respectively. Spoligotype analysis 
revealed that M. bovis strains from lechwe and cattle in Kafue basin clustered into a major node SB0120, where 
isolates outside the Kafue basin clustered into different nodes of SB0131 and SB0948. The comparatively higher 
variety of strains in cattle compared to lechwe elucidated by Mycobacterial Interspersed Repetitive Units-Variable 
Number Tandem Repeats analyses are consistent with cattle being the probable source of M. bovis in wild and 
domestic animals interphase area in Zambia. 

Conclusions: These results provide strong evidence of M. bovis strains transfer between cattle and lechwe, with the 
latter having developed into a sylvatic reservoir host. 
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Background 

Majority of the mycobacterial species that cause human 
and animal tuberculosis are grouped together as mem- 
bers of the Mycobacterium tuberculosis complex (MTC) 
[1,2]. This Mycobacterium tuberculosis complex includes 
very closely related species of mycobacteria among them: 
M. tuberculosis, M. africanum, M. microti, M. bovis, M. 
caprae, M. canettii and M. pinnipedii [1]. Although M. 
tuberculosis infection is the most common cause of 
human tuberculosis, part of other proportion of cases 
are due to M. bovis [3,4] . 

Zoonotic tuberculosis is caused mainly by M. bovis 
that has been shown to have a very wide host range 
[4-6]. The specie has been documented throughout the 
world with a similar impact in terms of disease 
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occurrence [5]. In Zambia, BTB is not homogenously 
distributed, however, high prevalence rates have been 
recorded within and around the Kafue basin where 
there is extensive overlap in terms of grazing land for 
both wild and domestic animals [7-9]. Additionally, the 
Kafue lechwe antelopes (Kobus leche Kafuensis) found 
in the Kafue basin have been described as feral reser- 
voirs of BTB in Zambia [10,11]. The disease has a his- 
torical presence in the Kafue basin that predates the 
identification of the area as a protected ecosystem and 
Ramsar Site no.530 [12]. Despite the continued reduc- 
tion in annual rainfall figures under the effects of glo- 
bal warming, the Kafue basin still remains as one of 
the few lucurstrine wetland ecosystems in Zambia and 
Africa, supporting a surging cattle population esti- 
mated at 300,000 animals, [13] at a carrying density of 
50 animals per square kilometre and approximately 
38,000 lechwe antelopes [12] on a 6,000 square kilo- 
metre wetland [14]. The area is characterised by high 
BTB with a herd level prevalence of around 50%, 
whereas a comparatively lower herd prevalence aver- 
aging 5.6% has been determined in areas outside the 
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basin [7,15]. Likewise, the corresponding Kafue lechwe 
antelopes have been shown to have a higher BTB 
prevalence rate [11,16], raising questions on a possible 
interspecies transmission of the disease between cattle 
and Kafue lechwe antelopes. This is however hampered 
by the lack of direct evidence to conclusively ascertain 
this assertion. 

Sequencing of the whole genome of the members of 
MTC [17-21] has shown a high level of sequence homo- 
geneity among the members (99.95%). Thus a careful 
and detailed comparative exploration into the individual 
genomes of the members of this complex was employed 
to mine out significant differences for diagnostic capabil- 
ities. Comparative genome analysis informed us that M. 
bovis has a smaller genome compared with M. tubercu- 
losis [1]. Furthermore, M. bovis has over time lost off 
some genes compared to M. tuberculosis. These genomic 
insertion-deletions are commonly referred to as Regions 
of difference (RD) and have been used in speciation of 
members of this complex as well as in explaining the evo- 
lution of the MTCs [1,17,19]. Spoligotyping diagnostic 
technique highlights intra species differences determined 
by the loss of spacers at a direct repeat region in MTCs, 
thereby creating a fingerprint typical of a particular specie 
[22]. Additionally, it is a more rapid and specific method 
of MTC speciation apart from being less laborious com- 
pared with biochemical, phenotypic and IS6110- restric- 
tion fragment length polymorphism (RFLP) analysis [22]. 
Spoligotyping results are very practical and reproducible 
across different laboratories internationally [22,23]. The 
technique has developed considerably and it is marked by 
a system of nomenclature and strain data capture and 
identification with a huge geographical and epidemio- 
logical relevance worldwide [24] . Most strains of M. bovis 
have one copy of IS6110 and spoligotyping is in general 
more discriminative when used with methods based on 
PCR amplification of the loci containing variable number 
tandem repeats (VNTRs) [25,26]. 

The target of this study was to molecularly characterize 
a population sample of M. bovis from cattle and kafue 
lechwe antelopes in Zambia to determine the genetic di- 
versity and relatedness of the isolates from domestic ani- 
mals and wildlife. 

Results and discussion 

Isolation and confirmation of M. bovis by MTCD-MPCR 

A total of 315 specimens from cattle and 75 from 
lechwe antelopes were analysed to initially determine the 
prevalence of MTC species. The samples were collected 
based on observations of gross pathological lesions upon 
examination. The observations included generalized lesions 
involving the lungs, pleural and mediastinal lymph nodes 
in both lechwe and cattle. These tuberculous lesions were 
observed by other workers who were investigating gross 



pathological distribution of tuberculous lesions in both cat- 
tle and lechwe [9,27,28]. From the samples analysed, 52 
MTC strains (Table 1) were obtained with 37 isolates from 
315 cattle and 15 from 75 lechwe. Following isolation M. 
bovis was confirmed by screening using the MTCD-MPCR 
analysis as previously observed [11,29]. The MTCD-MPCR 
technique is very useful in the differentiation of MTC at 
whatever level of diagnosis as it is simple and specific [29- 
31]. The presence of M. bovis in wildlife may translate into 
a perpetual focus of the disease [10] considering that BTB 
control in wild animals is a very difficult undertaking 
[32,33]. The only workable solution would be to intensify 
BTB testing in domestic animals so that reactors are 
removed. Furthermore traditional cattle herders must also 
be informed of the dangers of grazing their animals in 
areas where lechwe antelopes are present. 

Spoligotyping and Multiple locus variable number of 
tandem repeats analysis 

Spoligotyping of the 52 M. bovis isolates revealed their 
molecular clonality (Table 2). Two major spoligotype 
patterns (SB0120 and SB0131) were observed accounting 
for 36 isolates (69.2%), and 15 isolates (28.8%) cattle and 
lechwe respectively. Two isolates not identifiable with 
the two major clusters was given an SB0948 under the 
global spoligotype patterns diversity provided by the 
international data base on spoligotyping [24]. This only 
accounted for 4% of the observed isolates. M. bovis 
strains from both cattle and lechwe of the Kafue basin 
were found to share the same spoligotype (SB0120). This 
spoligotype was previously reported to be the major 
strain circulating in cattle around the Kafue basin, al- 
though by that time, no strains were determined from 
wild animals [34]. M. bovis strains isolated outside the 
Kafue basin were found to share a different cluster 
(SB0131). All the 52 isolates lacked spacers 3, 9, 16 and 
39 to 43, a characteristic feature that distinguishes M. 
bovis from M. tuberculosis [22]. Molecular studies have 
previously shown that clonality implies active transmis- 
sion of disease [35]. Thus our results in the Kafue basin 
where there was high clonality between isolates from 
lechwe and cattle, suggest an active transmission of M. 
bovis between the two animal species. The spoligotyping 
results showed that SB0131, clustered outside the Kafue 
basin suggesting that none of the lechwe antelopes har- 
boured M. bovis with such a spoligo pattern. 



Table 1 Results of the MTCD-MPCR of the isolated 
Mycobacterial isolates from cattle and lechwe 



Animal species 


Total sampled 


MTCD-MPCR M. bovis +ve 


Cattle 


315 


37 


Lechwe antelope 


75 


15 


Total 


390 


52 
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Table 2 Spoligotyping results by area of origin across animal species 


Animal species 


Area of origin 


Strain type 


Clade 


STI 


No of isolates 


Lechwe 


Blue lagoon National park 


SB0120 


BOVIS 1_BCG 


482 


4 




Lochinvar National park 


SB0120 


BOVIS 1_BCG 


482 


11 


Cattle 


Namwala (Southern Zambia) 


SB0120 


BOVIS 1_BCG 


482 


20 






SB0131 


BOVIS 1 


594 


4 




Lusaka (Central Zambia) 


SB0131 


BOVIS 1 


594 


1 




Chongwe (Central Zambia) 


SB0120 


BOVIS 1_BCG 


482 


1 






SB0131 


BOVIS 1 


594 


1 




Kabwe (Central Zambia) 


SB0131 


BOVIS 1 


594 


8 




Mumbwa (Central Zambia) 


SB0948 


No data 


New 


1 




Nampundwe (Central Zambia) 


SB0948 


No data 


New 


1 


Total 










52 



MLVA analysis identified 2 distinct clusters, which are 
corresponding to either SB0131 or SB0120 (Figure 1). 
Within the cluster possessing the spoligo pattern 
SB0131, there were 2 groups and 1 strain distinguishable 
one another by 2 loci, MIRU23 and QUB3232. The clus- 
ter possessing the Spoligotype SB0120 was differentiated 
in 6 groups. Nineteen strains, 14 from lechwe and 5 
from cattle, consisted of a major group and other 5 
groups or singletons were discriminated from the major 



group by a single MLVA locus difference each. The spo- 
ligotype SB0948 seems to be a descendant pattern of 
SB0120 lacking spacer 1 and also clustered in SB0120 
group with 1 locus (MIRU24) difference. The high isola- 
tion frequency of spoligotype SB0120 in cattle and 
lechwe conforms to the finding of previous works from 
Zambia [34] and may demonstrate dominance of this 
spoligotype despite the small sample size. The SB0120 
spoligotype has a considerable degree of geographical 
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Figure 1 Dendrogram based on the M. bows MLVA clusters distribution data. The major spoligotype SB0120 observed in cattle and lechwe 
are presented with varying MLVA types. The spoligotype SB0131 is also indicated with 3 MLVA types. 
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dispersion in Zambia and other African countries such 
as Algeria [36] and South Africa [37] . This spoligotype is 
also common in continental Europe [37-39]. Such a 
strain could have been introduced through livestock de- 
velopment schemes of trying to improve African beef 
local breeds. This finding has been highlighted by differ- 
ent workers [36,40,41]. Spoligotype SB0131 has been 
reported in this work for the first time in Zambia. This 
spoligotype has been documented in South Africa [42], 
which has a close livestock trade link with Zambia. 

To our knowledge, this is the first study conducting 
molecular characterization of M. bovis strains from 
cattle and lechwe in Zambia. However, routine studies 
have already previously revealed the presence of BTB 
in these animals [7,11,16,34]. In this study, results have 
shown genetic relatedness between M. bovis in cattle 
and lechwe antelopes of the Kafue basin. It is import- 
ant to note that the area around Lochnivar National 
park was used for ranching purposes by the early set- 
tlers from South Africa. This area was only gazetted as 
an animal sanctuary in 1972. It is the only place where 
Kafue lechwe antelopes antelopes are confined. The 
animals are semi-aquatic living in large groups near 
water and as such are confined to the area shown in 
Figure 2 (dotted circle). 

Characterization of M. bovis strains based on molecu- 
lar tools is important in understanding the epidemiology 
of BTB [25,36,41]. These results are significant in under- 
standing the transmission and dispersion of M. bovis 



strains within Zambia (Figure 2), given the high level of 
internal migration by the local people from the Southern 
part of the country to the Central and Northern regions 
of Zambia with their cattle. This type of internal migra- 
tion may lead to the dispersion of the SB0120 strain 
which in this present study was found confined to the 
Southern regions of Zambia. 

Conclusion 

The current study has described the possible source of 
M. bovis in wildlife and the transmission of limited 
strains of M. bovis between cattle and lechwe. The iden- 
tified maintenance and spread of M. bovis may become a 
dynamic and highly active process considering that 
human activities allow movement of cattle from one lo- 
cation to another. 

Methods 

Study area and sampling 

All the samples used in this study were collected in 
Zambia from 2006 to 2010. The bulk of the samples 
were obtained at abattoirs from cattle of the Kafue Basin 
and from other areas supplying cattle to the main abat- 
toirs of Lusaka, the capital city of Zambia (Figure 3). 
The animals were coming from the Southern province 
and Central areas of Zambia. Carcasses were followed 
along the slaughter line and screened for any visible tu- 
berculous lesions. All animals suspected to have tuber- 
culous lesions were examined further by collecting 





200 KM 



Figure 2 Spread of M. bow's strains in Zambia. Circle and dotted circle indicates the M. bovis from cattle and lechwe, respectively. White 
portion in circles indicates the M. bovis isolates with spoligotype SB0120, Black portion of circles indicates the M. bovis isolates with spoligotype 
SB0131, and the shaded circle indicates the M. bovis with spoligotype SB0948. Size of circle indicates the strain number. 
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Figure 3 Study area in Zambia. A: Blue lagoon National park, B: Lochnivar National park, C: Namwala, D: Lusaka, E: Chongwe, F: Kabwe, G: 
Mumbwa, H: Nampundwe. A, B, C, G and H are an interphase area consisting of a wildlife sanctuary and a traditional cattle grazing area; while 
Area D, E and F are a farming area consisting of commercial and subsistence cattle farms. 



samples for laboratory analysis. The wildlife samples 
were obtained from Kafue lechwe antelopes, which were 
sampled over the same period as cattle through a special 
research license from the Zambia Wildlife Authority. 
The sample size depended on the number of animals 
allocated for scientific studies with the selection being 
done at random by gun stunning. The animals were col- 
lected from Blue lagoon national park (Figure 3: area A) 
and Lochnivar national park (Figure 3: area B). The ani- 
mals were examined for gross lesions according to the 
standard post mortem examination procedures [43]. The 
organs and tissues with suspected TB lesions were col- 
lected together with accompanying demographical data 
of area of origin, type of organ or tissue as well as the 
type of gross pathological postmortem disposition. These 
specimens were placed in sterile self-zipping histopatho- 
logical bags, placed into a cooler box with ice packs be- 
fore transportation to the laboratory for analysis. 

Cultivation of Mycobacterial and DNA extraction 

The collected samples of TB suspect were analyzed and 
cultured for growth as previously described [16]. Briefly, 
the tissues were trimmed of fat and then a 500 mg sam- 
ple was minced with sterile scissors and homogenized in 
a sterilized glass homogenizer. A milliliter of phosphate 
buffer (pH 6.8) was added and thoroughly mixed after 
which 1 mL of 5% sodium hydroxide was added. After 
incubation for 15 min at room temperature, 10 mL of 
phosphate buffer was added and then centrifuged at 



1500 xg for 20 min. The pellet was collected and then 
resuspended in a final volume of 0.5 mL of phosphate 
buffer which was then used for inoculation onto 2% 
Ogawa medium. Bacterial growth was then monitored 
up to 8 weeks at 37°C. The resulting cultures were tenta- 
tively identified as probable Mycobacterium tuberculosis- 
complex by their slow growth and colony morphology. 
Purity and acid-fastness of the colonies were checked by 
Ziehl Neelsen staining. DNA was extracted from Myco- 
bacterial colonies using DNAzol reagent (Invitrogen, 
Carlsbad, CA, USA) and mechanical disruption as previ- 
ously described [44] according to the manufacturers 
instructions and dissolved in 50 uL TE buffer consisting 
of 10 mM Tris/HCL (pH 8.0) and 1 mM EDTA. 

MTC discrimination by multiplex PCR 

The Mycobacterium tuberculosis complex-discriminating 
multiplex PCR (MTCD-MPCR) targeting genetic regions 
cfp32 (a specific gene for MTC), RD9 (region of difference 
9; seen only in M. tuberculosis and M. canettii) and RD12 
(region of difference 12; deleted in M. bovis, M. caprae 
and M. canettii) was used for species differentiation 
according to the previously publication [27] . The reaction 
mixture consisted of 7.4 uL H 2 0, 2 uL 10 x Taq buffer, 2 
uL dNTPs (2.5 mM each), 0.2 uL Taq (Takara Bio Inc, 
Shiga, Japan), 1 uL of DNA sample, 2.2 uL of 10 uM cfp32 
primers, 0.7 uL of 5 uM RD9 primers and 0.8 uL of 5 uM 
RD 12 primers. The PCR was performed using the follow- 
ing program: denaturation for 1 min at 98°C followed by 
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35 cycles of 5 sec at 98°C, 20 sec at 58°C and 1 min at 
68°C with final elongation at 72°C for 5 min in a thermal 
cycler (iCycler, Bio-Rad Laboratories Inc., Hercules, CA). 
The PCR products were separated by electrophoresis in a 
2% agarose gel in TAE buffer, and then visualized after 
staining with ethidium bromide. 

Spoligotyping of M. bovis isolates using micro-spoligoarrays 

Spoligotyping was performed according to the procedure 
by Kamerbeek and co-workers with slight modifications 
[22]. Forty- three spacer-sequence probes were covalently 
bound to the membrane (Pall Co., NY, USA). The pri- 
mers used were DRa (GGTTTTGGGTCTGACGAC) 
and DRb (CCGAGAGGGGACGGAAAC). A hot start PCR 
was done by mixing 1 uL each of the 10 uM primer, 7.3 uL 
H 2 0, 3 uL of 5 x colorless Go Taq buffer (Promega™, 
Fitchburg, WI), 1.5 uL of PCR DIG Labeling Mix (Roche), 
0.2 uL of Go Taq DNA Polymerase (5 units/ uL; Pro- 
mega) and 1 uL of DNA sample in 15 uL total reaction 
mix per tube. The PCR reaction was initiated by de- 
naturation at a temperature of 98°C for 1 min, followed 
by 40 cycles of 98°C for 5 sec, 55°C for 10 sec and 72°C for 
30 sec with a final elongation step at 72°C for 1 min in a 
thermal cycler. The 500 times diluted PCR product in 
hybridization buffer was heat denatured at 95°C for 5 min- 
utes and immediately cooled on ice to leave the DNA sin- 
gle stranded. Hybridization was performed by placing the 
nitrocellulose membranes (Pall Co., NY, USA) for incuba- 
tion at 60°C for 1 hour. After which, the membrane was 
washed in TBST-E (0.1% Tween-20 and ImM EDTA-2Na 
in TBS) at 60°C for 1 min and thenlO min, finally 1 min. 
The membrane was then dried at room temperature. DIG 
on the nitrocellulose membranes were reacted with a 1000 
times diluted Ant-Digoxigenin-POD [poly], Fab fragments 
(Roche), with TBSTE-E at room temperature for 30 min. 
Then the membranes were sequentially washed in TBST-E 
at room temperature for 1 min, 10 min and 1 min. Then, 
POD on the membranes was detected by TMB solution 
(TMB Peroxidase Substrate Kit, Vector Labs Inc™, Burlin- 
game, CA) according to the manufacturer s protocol. 

Multiple locus variable number of tandem repeats (MLVA) 
assay of M. bovis isolates 

The isolates were further genotyped by PCR amplifica- 
tion using primers targeting 26 VNTR loci [45]. Two 
different PCR reaction mixtures were conducted accor- 
ding to the loci. Locus MIRU10, 16, 24, 26, 27, 39, 40, 
ETR-B, F, VNTR-424, 3690 were done under betaine 1.0 
M, whilst locus MIRU2, 4, 20, 23, 31, ETR-A, C, QUB- 
11a, lib, 26, 3232, 3336, VNTR-1955, 2401, 4156 were 
performed under GCII buffer (Takara). The PCR reac- 
tion mixture under 1.0 M betaine buffer was conducted 
in a mixture consisted of 6.3 uL H 2 0, 0.4 uL of dNTP 



(10 mM each), 3.0 uL of 5x colorless Go Taq buffer, 3.0 
uL of betaine (5 M), 0.6 uL of Primer 1 and 0.6 uL of 
Primer 2, 0.1 uL of Go Taq (5 units/ uL) DNA Polymer- 
ase, and 1 uL of DNA sample. An initial denaturation 
step of 95°C for 5 min was followed by 32 cycles at 95°C 
for 15 sec, 58°C for 20 sec and 72°C for 1 min with a 
final elongation step at 72°C for 1 min in a thermalcy- 
cler. The PCR reaction mixture under the GCII buffer 
consisted of 4.8 uL H 2 0, 0.4 uL of dNTP (10 mM each), 
7.5 uL of 2 x GCII buffer, 0.6 uL of primers, 0.1 uL of 
Go Taq (5 units/ uL) DNA Polymerase (Promega), and 1 
uL of DNA sample. An initial denaturation step at 95°C 
for 5 min was followed by 32 cycles of 95°C for 15 sec, 
50°C for 20 sec and 72°C for 45 sec with a final elong- 
ation step at 72°C for 1 min in a thermalcycler. All the 
samples were electrophoresed on a 2% agarose gel to 
identify the repeat numbers. 
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